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Abstract. Dynamics of molecular motor dimers, consisting of rigidly bound particles 
that move along two parallel lattices and interact with underlying molecular tracks, is 
investigated theoretically by analyzing discrete-state stochastic continuous-time burnt- 
bridge models. In these models the motion of molecular motors is viewed as a random 
walk along the lattices with periodically distributed weak links (bridges). When the 
particle crosses the weak link it can be destroyed with a probability p, driving the 
molecular motor motion in one direction. Dynamic properties and effective generated 
forces of dimer molecular motors are calculated exactly as a function of a concentration 
of bridges c and burning probability p and compared with properties of the monomer 
motors. It is found that the ratio of the velocities of the dimer and the monomer 
can never exceed 2, while the dispersions of the dimer and the monomer are not very 
different. The relative effective generated force of the dimer (as compared to the 
monomer) also cannot be larger than 2 for most sets of parameters. However, a very 
large force can be produced by the dimer in the special case of c = 1/2 for non-zero 
shift between the lattices. Our calculations do not show the significant increase in the 
force generated by collagenase motor proteins in real biological systems as predicted 
by previous computational studies. The observed behavior of dimer molecular motors 
is discussed by considering in detail the particle dynamics near burnt bridges. 
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1. Introduction 

Theoretical and experimental investigations of motor proteins, or molecular motors, 
have recently become a subject of considerable attention because of importance of these 
molecules for fundamental understanding of non-equilibrium processes in chemistry and 
biology [H El El Hj. Motor proteins are active enzymatic molecules that consume a 
chemical energy and transform it into mechanical work, while moving at the same time 
along linear molecular tracks. They play a critical role in many biological processes 
[DOE]- In contrast to conventional molecular motors powered by the hydrolysis of 
ATP or related compounds, a recently discovered molecular motor called collagenase 
fuels its motion along collagen fibrils by using a significantly different mechanism that 
utilizes asymmetric collagen proteolysis [51 [6] . The biased diffusion of this motor protein 
along the molecular track results from the fact that the motor protein molecule after 
cutting the bond at the specific site of the collagen fiber is always found on one side of 
the reaction site, and it is unable to cross the broken link after the cleavage. 

Analysis of collagenase transport suggests that a successful description of this motor 
protein's dynamics can be obtained by utilizing so-called "burnt-bridge models" (BBM) 
PS EJ El El [in] ■ According to BBM, the molecular motor is described as a random 
walker hopping along the one-dimensional lattice composed of two kinds of links: strong 
and weak. The walker does not affect the strong links when crossing them, while the 
weak links ("bridges") might be destroyed ("burnt") with a probability < p < 1 
when crossed by the walker; the walker cannot pass over the destroyed links again. In 
Refs. [9l [10] we developed a new analytical approach that allowed us to obtain exact 
expressions for motor protein velocity V(c,p) and dispersion D(c,p) for different sets of 
parameters and for different burnt-bridge models, and it was shown that the analytical 
predictions are in full agreement with Monte-Carlo computer simulations. 

Theoretical and experimental investigations of collagenase dynamics [S, E] conclude 
that this motor protein is rather weak: at standard conditions it exerts a stall force of 
Fs — 0.02 pN, which is 100 times smaller than other conventional motor proteins [4], 
and thermodynamic efficiency of this motor protein is between 6 and 15%, although it 
has a significantly larger processivity than conventional motor proteins [21 E]. Meantime, 
experimental studies of several collagenases capable of breaking collagen fibers [6lfTTl[T2] 
indicate that in cells these enzymes are probably clustered together. Computer 
simulations of the motion of rigidly bound dimers [6] of collagenase molecules produced 
stall forces of the order of 5 pN, 10-100 times larger than the single motor protein, 
and it was suggested that cell surface-anchored clusters of collagenases might play an 
important role in cell motility. Coupling of motor proteins can produce very efficient 
molecular motors and it can lead to interesting dynamic phenomena [131 EH] • However, 
the motion of dimer particles in the context of BBM has not been studied yet. The 
goal of the present work is to investigate theoretically the dynamics of dimer molecular 
motors in BBM using previously developed exact analytical approaches 0, [JO] • 



Transport of Molecular Motor Dimers in Burnt-Bridge Models 
2. Model 



3 



We consider a molecular motor dimer as two particles connected by a rigid bond that 
move along two parallel infinite one-dimensional lattices (a bond between the particles 
is perpendicular to the linear molecular tracks), with transition rates in both directions 
and the lattice spacing assumed to be equal to one, as shown in Fig. 1. The links 
connecting consecutive binding sites on the lattices fall in two categories. While the 
random walkers have no effect on the strong links, each of them destroys the weak ones 
with the probability p when passing over them. After the weak link is broken, the 
particle is always assumed to be to the right of it, and all weak links are intact at t = 0. 
Therefore, the dimer cannot be trapped between two burnt bridges, and it continuously 
moves to the right. 
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Figure 1. Schematic view of the motion of molecular motor dimer in the burnt-bridge 
model. Strong links are represented by thick solid lines, while thin solid lines depict the 
weak links which can be destroyed. Dotted lines correspond to already burnt bridges. 
The dimer can jump with equal transition rates one step to the right or to the left 
provided that the link is not already destroyed. The bridges on the upper and lower 
tracks arc shifted by I lattice spacings. 

The bridge can be burnt according to one of the two scenarios. In the first 
mechanism, the bridge burns only when the random walker crosses it from left to right, 
however it always remains intact if crossed from right to left. We term this model a 
forward BBM. According to the second model, the weak link gets broken if crossed in 
either direction, although the particle is always found to the right of the broken bond. 
This is a forward-backward BBM. The physics of the driven motion is similar for both 
mechanisms, and both models are even identical for p — 1, however for p < 1 the walker's 
dynamics in two cases are different [9] . In this paper we will consider only forward BBM 
because theoretical analysis is simpler in this case. However, our approach can be easily 
extended to more complex forward-backward BBM. 

The dynamics of the dimer in BBM depends on two parameters: the burning 
probability p and the concentration of bridges c. The properties of the system are 
also strongly influenced by the distribution of bridges, periodic or random [8J. In the 
case of the periodic distribution, which will be analyzed in detail in the present paper, 
the bridges on each of two parallel one-dimensional lattices are positioned equidistantly, 
with a constant distance N = 1/c lattice spacings separating them. This distribution 
is more suited for the analysis of the transport of collagenases since the sites that can 
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be destroyed by enzymes on collagen fibers are found to be at equal distances from 
each other [5]. All the models presented below will be studied in the continuous-time 
description that better corresponds to the real biological systems with motor proteins 
[9]. We also assume that the bridges on two parallel tracks are shifted by I lattice 
spacings (0 < / < N — 1): see Fig. 1. 



3. Dynamics of Molecular Motor Dimers 

3.1. BBM withp= 1 
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Figure 2. Schematic picture for the transport of the single particle (equivalent to the 
dimer in BBM with p = 1) along the effective single linear track. 



Let us consider first a special case of BBM with p = 1, i.e., when a bridge burns 
every time when crossed by the walker. It can be easily seen that in this case the motion 
of a molecular motor dimer along two parallel tracks is equivalent to the motion of a 
single random walker on a single track but with increased concentration of bridges [6], 
as shown in Fig. 2. The distribution of bridges on this effective track is again periodic 
but generally non-uniform, with two weak links in each period: between the sites I — 1 
and /, and between the sites N — 1 and N: see Fig. 2. The dynamics of this system can 
be mapped into a random single-particle hopping model on periodic lattices that was 
solved by Derrida [U [15], in which the explicit expressions for the dynamic properties 
are obtained. We define the forward and backward transition rates for the random 
walker at site j (j = 0, 1, • • • , N — 1) as Uj and Wj, correspondingly. The velocity and 
the diffusion constant are then given in terms of these transition rates [H IT5] . 
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The transition rates for the particle shown in Fig. 2 are Uj = Wj = 1 for all j except 
w = wi = (the burnt bridges cannot be crossed again). Then the explicit expressions 
for (J3]) and (j4j) are the following, 

for j = 0, 1, 1; 








for j = 1,1 + 1,-- -, JV- 1; ^ 



and 



_/ J + l, for j = 0, 1, — 1; 



j-Z + 1, for 7 = Z,Z + 1,---,JV-1 

Substituting these results into Eq.([T]) leads to the expression for the dimer's velocity, 

2JV 

N) = i(i + i) + (N-i)(N-i + iy (7) 

Calculation of the dispersion D according to (J2]) is more complex, but it can be 
performed with the help of the following useful identities: 

+ i + for < i < I - j - 2; 

'-,-,-,+.! ^ iV-(j + z + l), for/-j-l<i<JV-j-2; (8) 
Z + JV - (j + i + 1), for N - j - 1 < i < N - 1, 

for j < Z — 1, while for j > Z we have 

f JV - (j + i + 1), for < % < N - j - 2; 
r,- +i+ i = | / + N - (j + i + 1), for JV - j - 1 < i < N + I - j - 2; (9) 
[ 2JV — (j + i + 1), for JV + Z - j - 1 < z < N - 1. 

As a result, the dispersion for the dimer is given by the following formula, 

9/V 2 r 

= 3(21' + AT - 2IAT + Af)» ^ " " N ~ 2W(2 + ™ + ^ > 
+2Z 2 (2 + 3JV + 3JV 2 ) + JV(1 + JV)(1 + JV + JV 2 )J . (10) 

The validity of Eqs. (|7j) and (flOj) is confirmed by the fact that for Z = and Z = JV 
they reproduce the BBM results obtained for the transport of monomers on the lattice 
with the period JV j9], 

y= J_ fl= 2(Jf + JV+l) . 

JV+1' 3(/V + l) 2 ' 1 J 

whereas for N = 21 we recover the results for the monomer's dynamics on the lattice 
with the period Z, as expected. Also, Eqs. (j7]) and ffTO]) are invariant under the change 
/ — > JV — Z, reflecting the symmetry of the system. 

For the purposes of comparison of the dynamic properties of the dimer with those 
of the monomer in BBM with p — 1, we plot their ratios in Figs. 3 and 4 as functions 
of the concentration of bridges c = 1/JV and the shift I between the parallel tracks. The 
range of Z values (for constant c) is < I < JV = 1/c, whereas that of c values (for 
constant I) is < c < l/Z. Hence it is convenient to express the dynamic properties of 
the dimer as functions of dimensionless parameter Ic (with either fixed Z or c), so that 
< lc < 1 and the curves with different I (or c) values could be plotted on the same 
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Figure 3. (a) Ratio of dimer's and monomer's velocities as a function of the 
concentration of bridges c for the fixed value of the shift I. (b) Ratio of dimer's and 
monomer's velocities as a function of the shift I for the fixed concentration of bridges 
c. 
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Figure 4. (a) Ratio of dimer's and monomer's dispersions as a function of the 
concentration of bridges for the fixed value of the shift I. (b) Ratio of dimer's and 
monomer's dispersions as a function of the shift I for the fixed concentration of bridges 
c. 



interval. The physical meaning of the parameter Ic is a dimensionless concentration of 
bridges (for fixed I), or a dimensionless shift distance (for fixed c). The invariance of 
dynamic properties © and (jTOl) under replacement I — > N — I implies that when plotted 
as functions of I (for fixed c), their graphs are symmetric with respect to the middle of 
the Ic interval, i.e. lc=\: see Figs 3b and 4b. 

As shown in Fig. 3a, the relative velocity of the dimer has a non-monotonic behavior 
as a function of the concentration of weak links, and V^/Vi reaches a maximum at 
c* = 1 at which V 2 (l, c*)/Vi(c*) = 1 + - J[ < 2. For very large shifts, I > 1, 
one can easily see that c* — > h an d V^il, c*)/Vi(c*) ~ 2 — y. Similar behavior is observed 
for the relative velocity of the dimer as a function of the shift distance Z, as illustrated in 
Fig. 3b. As discussed above, because of the symmetry the maximum ofV^/Vi is observed 
at I* = y = ^, at which V%{1* c]/V\{c) = ? an d f° r ver Y l° w concentration of 

bridges (N — > oo) the dimer molecule moves twice as fast as the monomer molecular 
motor at these conditions. Thus, we conclude that at all conditions V 2 (l, c)/Vi(c) < 2. 
It is also worth mentioning that the dimer's velocity itself V-^/, c) reaches the maximum 
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at c = -jy (for fixed I) and I = I* = y = y (for fixed c), as can be deduced from Eq. 

([7]). The corresponding values of V 2 are V^/, c) = 2 (y/2-i)i+i an< ^ ^(^' c ) = W+2- 

The dispersion of the dimer shows a more complex behavior with two maxima and 
one minimum, as presented in Fig. 4. Our calculations show that generally D 2 /D\ lies 
between 0.92 and 1.31. It is interesting to note that for some range of parameters (not 
very large shifts and not very small concentration of bridges) the dimer fluctuates less 
than the monomer, leading to D 2 /Dx < 1. 

Simultaneous knowledge of velocities and dispersions of the dimer allows us to 
evaluate the effective force exerted by the molecular motor [9] . The motion of the dimer 
in the system with burnt bridges can be viewed as an effective biased random walk in 
the system without bridges. Then its dynamic properties can also be written as 

V = a{u eff - w eff ), D = ^a 2 {u eff + w eff ), (12) 

with u e ff and w e // being the effective transition rates of the biased random walk [9]. 
The effective force generated by the dimer is given by [U [16] 



F= ksT ln Ue iL = ksT ln 2D + Va^ 
a w eff a 2D-Va 
where a is the size of the discrete steps made by the walker, and it corresponds to 
the distance between the consecutive lattice sites. Parameters V and D are expressed 
in units of rajs and m 2 /s correspondingly, and they are related to the dimensionless 
dynamic properties V and D [Eqs. (J7|) and ffTUj) ] in the following way: 

V = av V, D = a 2 v D, (14) 

where Vq is the intrinsic transition rate of the molecular motor. Then 

with functions V(l,N) and D(l,N) given by Eqs. ([7]) and fflOl respectively. 

To investigate the efficiency of the dimer molecular motor we compare the force 
generated by the dimer particle in BBM with the force exerted by the monomer at the 
same conditions, and results are presented in Fig. 5. Note that F 2 (l,c) is given by Eq. 
f fT5|) . while the corresponding force for the monomer is equal to i*i(c) = F 2 (l = 0, c). As 
shown in Fig. 5a, for / > 2 the ratio of effective forces is less than 2 for all concentrations 
of bridges, although the maximum of the ratio might slightly exceed 2 for small / (~ 2.04 
for 1 = 2). The situation is drastically different for Z = 1 where F2/F1 — > 00 as c — ► 1/2 
(since lim F 2 = 00 for 1 = 1). For all curves in Fig. 5a, the maximum is reached 

c-i-l/2 

very close to Ic = 0.5, and it approaches Ic = 0.5 as / — > 00. For fixed concentration of 
bridges (see Fig. 5b), symmetry requires that F 2 (N — l,c) = F 2 (l, c), and F 2 /Fi reaches 
its maximum at I* = y = 5-. It is clear from Fig. 5b that except for the case of c = 1/2, 
where F 2 jF\ ^00 as I — > 1, the maximum of the ratio F 2 jF\ is never larger than 2. It 
approaches 2 (from below) as iV — >• 00 (c — >• 0). 
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Figure 5. (a) Ratio of dimer's and monomer's effective forces as a function of the 
concentration of bridges c for the fixed value of the shift I. Note that for I — 1 there 
is no branch for Ic > 1/2 because there are no possible values of c between 1 and 1/2 
(N is integer), (b) Ratio of dimer's and monomer's effective forces as a function of the 
shift I for the fixed concentration of bridges c. 



3.2. BBM with p < 1: a special case 

Although the general situation when the burning probability is p < 1 requires a separate 
consideration that leads to extensive calculations (see next section), there is a special 
case of BBM for which the dynamics of the dimer molecular motor can be expressed in 
terms of the properties of the effective monomer in BBM, for which analytical expressions 
are already known [9], [10]. When bridges on the parallel tracks are shifted by the distance 
I = N/2, i.e., half of the period, the motion of the dimer is identical to hopping of the 
effective particle on the single lattice with the distance of N/2 between neighboring 
weak links. Then, using results from Ref . [HI [10] , the relations between the dynamic 
properties of the dimer and the corresponding monomer are the following for I = N/2, 

V 2 (c,p) = V l (2c,p), D 2 (c,p) = D l (2c,p). (16) 

To simplify calculations we only consider forward BBM, in which bridges burn 
only when crossed from left to right, but remain unaffected if crossed in the opposite 
direction. For this model, analytical expressions for V\(c,p) and D\(c,p) are already 
obtained [9J [10]. The velocity of the monomer is 

V l {c,p)= 2P ° (17) 

lp 2 (l - c) 2 + Ape 



while the dispersion has more complex form [10], 
2 



-2A + 2 - R - 2R { — — — — — ^[iV + 1 + (N- l)x] - A * A 



(iV _2)(iv 2 -n 2 



24 u 1 

where parameters A and A are given by 




6(1 - a;) 



Ri - -A (18) 



= N{1 - x)[apx + Rq(1 -x)(-l+p + x)] 

p[l + (N- l)x](-l +p + x 2 ) ' 1 ' 
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Figure 6. (a). Ratio of dimer's and monomer's velocities as a function of concentration 
of bridges for I = l/(2c). (b) Ratio of dimer's and monomer's dispersions as a function 
of concentration of bridges for I = l/(2c). 



with 



and 



A = (l-p)[ap + R (l-x) 2 } 

p(-l+p + x 2 ) { ' 

a= r + {l-x)R , (21) 

r = ^(l -x)-Rl. (22) 

Thus Eq. ( I18p is expressed in terms of the parameters Rq and x which are also functions 
of N = 1/c and p, 

Ro = = = = = , (23) 

\Jp 2 (N - l) 2 + ApN 09 2 (1 - c) 2 + Ape 



x = 1 + -p(N — 1) — -yp 2 (N — l) 2 + ApN. (24) 

Therefore, Eq. ( [181) combined with Eqs. ( fl9l) - f l24l provides the explicit expression for 
Di (c, p) as a function of the concentration of bridges and the bridge burning probability. 

In this case the dynamic properties of dimer and monomer molecular motors can 
be compared for the entire range of burning probabilities, < p < 1, as presented in 
Fig. 6. Since the shift / assumes only integer values (I = 1, 2, • ■ ■), the relation le = 1/2 
implies that we only need to consider the set of discrete c values, c = 1/2, 1/4, 1/6,---, 
i.e., only the interval < c < 1/2 is relevant. As seen from Fig. 6a, decreasing the 
burning probability for any fixed concentration of bridges lowers the relative velocity 
of the dimer. Increasing the number of weak links (for fixed p) also lowers the relative 
velocity of the dimer, and the ratio of velocities lies within the interval, a/2 < V2/V1 < 2. 
This result can be explained using the asymptotic expressions for the velocity obtained 
earlier [TO] . 

K l( c,p)^{ * 1°'^/ (25) 
I Jpc, for pCc. 
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Then for c — > the ratio V2/V1 approaches 2 for any finite value of the burning 
probability, while for p <C 1 (for finite c) we have V2/V1 — y/2. As shown in Fig. 6b, the 



dispersion of the dimer is similar to the dispersion of the monomer, 0.92 < D 2 jD\ < 1 for 
< c < 1/2. Lowering the burning probability increases the relative diffusion constant 
of the dimer particle, while the effect of the bridge concentration is non-monotonic. 
Because in the limiting cases the dispersion is independent of the concentration of 
bridges [9] [Di(c,p) ~ 2/3 for c <C p and Di(c,p) ~ 1/2 for p < c], the ratio of 
dispersion is D 2 /D 1 — > 1 in these cases. 

The effective forces exerted by dimer and monomer molecular motors for p < 1 
and Ic = 1/2 can be obtained using the approach described in the previous section for 
the case of p — 1, and the results are plotted in Fig. 7. The force exerted by a single 
particle Fi(c,p) is found from Eq. (jl5l) with D = Di(c,p) and V = Vi(c,p) given by 
Eqs. (1181) and ( jTTl) . Similarly, the force generated by a dimer F 2 (c, p) is found from Eq. 
(|T5|) with D = D 2 (c,p) and V = V 2 (c,p), and thus F 2 (c,p) = Fi(2c, p). Making use of 
asymptotic relations for the velocity and dispersion [10] yields the following expressions 
for Fi in the limiting cases, 



It follows that F 2 jF\ — > 2 if c — > for any finite value of p, and F 2 /F\ — > \[2 for 
small burning probabilities (with fixed value of c). For p = 1, the ratio F 2 jF\ strongly 
increases as c — > 1/2. For p < 1, each curve shows a non-monotonic behavior with 
a minimum. The position of the minimum moves to the right with decreasing p, and 
it moves beyond the interval < c < 1/2 for p < 0.2. The minimal value of F 2 jF\ 
decreases with decreasing burning probability. The ratio F 2 jF\ exceeds 2 within the 
interval 0<c<l/2ifp> 0.74, this occurs in the vicinity of c = 1/2. In particular, for 
p close to 1 the ratio F 2 jF\ at c = 1/2 is substantially larger than 2 [e.g. at c = 1/2, 
F 2 /Fi w 2.52 for p = 0.9, F 2 /F 1 ta 3.94 for p = 0.99, F 2 /F 1 w 5.44 for p = 0.999]. Our 
analysis indicates that the relative force exerted by the dimer could substantially exceed 
2 only at c = 1/2 if the burning probability p is close to 1. 

Given that in the p = 1 case with fixed c the ratios V 2 /Vi and F 2 jF\ reach their 
maximal values at / = N/2 (as can be seen from Figs. 3b and 5b), it is reasonable to 
suggest that the same holds for p < 1. Thus the special case of Ic = \ considered in the 
present section provides the maximal values for ratios V 2 jV\ and F 2 jF\ for any given c, 
i.e. for I values other than I = 77- these ratios would be smaller. 

3. 3. BBM with p < 1 : general case 

Now let us consider the general case of the transport of the dimer molecular motors 
along two parallel tracks shifted by / lattice spacings with / = 0,l,2,---,iV — 1 and 
with the burning probability p < 1 (see Fig. 1). Again, we analyze the continuous-time 
forward BBM when the bridges are destroyed only when crossed from left to right. For 
burning probability p < 1 the Derrida's method is not applicable since we do not know 




for c<p; 
for p <^ c. 



(26) 
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Figure 7. Ratio of the force exerted by dimer and that exerted by monomer as a 
function of concentration of bridges c with shift I = 1/ (2c) . 

how the transition rates Uj and Wj are expressed in terms of p [9]. Instead, we will apply 
a method developed in Ref. [9] to investigate BBM with p < 1 for the monomer random 
walkers. It should be pointed out, however, that this approach allows us to calculate 
only velocities. To determine the diffusion constants one should utilize a more complex 
approach outlined in Ref. [ID] , but in this case the calculations become very involved 
and we will not discuss them here. 




Figure 8. Reduced kinetic scheme for the motion of the dimer in continuous-time 
forward BBM with bridges on the parallel tracks shifted by / lattice spacings. Only 
transition rates not equal to one are shown. Two origins correspond to the right end 
of the last burnt bridge at each track. 

Let us define Rj{t) [Uj(t)] as the probability that a random walker is located j 
sites apart from the last burnt bridge at time t given that the last burnt bridge lies 
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on the upper (lower) of the two tracks (see Fig. 1). The probabilities Rj(t) and 
Uj(t) are conveniently described in a reduced kinetic scheme presented in Fig. 8. The 
description of the dimer's dynamics in terms of the functions Rj(t) and Uj(t) corresponds 
to considering the system in the moving coordinate frame with the last burnt bridge 
always at the origin of the upper or lower linear track. The kinetic scheme in Fig. 8 
describes correctly the dynamics for I = 1, 2, • • • , N — 1. However, the case of zero shift, 
I = 0, requires a slightly different analysis, which will be given below. 

The time evolution of the system is determined by a set of Master equations. For 
the sites on the upper track the equations are 

kN ^ 1 ^ ' = RkN+i-l{t) + RkN+i+l{t) — 2RkN+i{t), (27) 

for k = 0, 1, 2, • • • and 1 < i < N — 1 (but with i ^ I). For the sites with % = I we have 

dRkN+i(t) 
dt " 

while for i = and k > 1 
dR kN (t) _ 
dt 

For the origin site the dynamics is more complex, 
dR Q {t) 



[1 — p)RkN+i-\{t) + RkN+l+i{t) — 2R kN+ i(t); (2£ 



1 -p)R kN ^(t) + R kN+1 (t) - 2R kN (t). (29) 



dt 



P 



X' 



RkN~i{t) + UkN-i-i(t) 
k=i k=i 



+ R 1 (t)-R (t). (30) 



The Master equations for Uj(t) are obtained from ( !27l) - ( J3D1 by changing R(t) — > U(t) 
and I —*■ N — I. Specifically, at the origin site of the lower track we obtain 

dU (t) 



dt 



V 



CO 



UkN~i(t) + 51 RkN+i-i(t) 

k=l k=0 



+ U 1 (t)-U (t). (31) 



The probabilities Rj(t) and Uj(t) should also satisfy a normalization condition, 

oo N-l 

E E [RkN+M + u kN+i {t)} = i. (32) 



fc=0 i=0 

In the stationary-state limit Eqs. (127)) - (1311 yield 



R i = R -ipE 1 , E x = E-^iV-i + E u kN-i-i, for < z < Z; (33) 

fc=i fc=i 

^AT+i+i = (i+l)i2 feJV +i-i(l-p)-RfeJV+i-i, for < i < N-l; k > 0; (34) 

-Rfc^+i = (i + l)i2fejv - «(1 - p)RkN-i, for < z < Z; fc > 1. (35) 

Similar results can be obtained for the functions Uj from Eqs. (1331) - (135)) by changing 
R —> U, I N — I, and Si — > E 2 , 

CO CO 

U i = U -ipT l2 , S 2 = E ^fev-i + E RkN+i-i, for < z < A^-Z; (36) 

fc=l fc=0 

U kN+N -i+i = (i+l)U kN+N -i-i(l-p)U k N+N-i-i, for < z < Z; A; > 0; (37) 
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U kN+i = (i + l)U kN - z(l - p)U kN ^, for < % < N - I; k > 1. (38) 



The systems of Eqs. ( 1331) - (1351) and Eqs. (1361) - (1381) are coupled via the auxiliary 
functions £i and E 2 . 

Following the approach of Ref. [9], we seek the solutions of Eqs. ( |33|) and ( |35l) in 
the form 

R kN+i = R e ak -iB(k), for < i < I, (39) 

where unknown parameters a and B{k) are to be determined. Expressing functions 
RkN+i and RkN+i-i with the help of Eq. fl39|) and substituting them into Eq. (134"|) . one 
finds: 

R kN+ i +i = R e ak -lB(k)+i [ P R e ak + B(k){-1 + p - pi}] for < i < N -1.(40) 

It is important to note that utilizing Eq. (|39|) to express R)~n+i-i implies that / — 1 > 0, 
i.e. / > 1. Next, we use the fact that Eq. (1401) with i = N — I should coincide with Eq. 
(1391) for i = 0, and the unknown function B(k) is eliminated, leading to 

B(k) = R e ak {l - e a + p(N - 1)}/{1 + (N - l)[l - p + pi]}. (41) 

Hence the system of Eqs. ( 1331 - ( |35l) is transformed into Eqs. ( 1391) and ( T^Ol) with the 
function B(k) given by Eq. (j4T|) . Similarly, the solutions of the system of Eqs. ( 1361) - 
(|38|) are sought in the following form, 

U kN+i = U e ak -iB(k), for <i < N -I. (42) 

Then we obtain 

U kN+N -i +i = U e ak -(N-l)B(k)+i [pU e ak + B(k){-1 +p-p(N -I)}], for < % < /(43) 

with 

B(k) = U e ak {l - e a + pl}/{N - I + l[l - p + p(N - I)}}. (44) 

To find the unknown parameter a, we utilize the recurrent formula (1351) with the 
functions R k N+i and R k N expressed via Eq. (1391) and the function R k N-i found from 
Eq. m, 

Roe ak - iB{k) = {i + l)R e ak - i(l - p) {Roe 11 ^ - IBik - 1) 

+[N-l-l] (pR a e a{k - l) + B(k - 1)[-1 + p - p/]) } . (45) 

Expressing -B(fc) from Eq. (T4T1) and introducing the parameter x = e a , one can reduce 
Eq. (|35D into 

a; 2 - [2 + 2p(A^- 1) +p 2 (/ - 1)(JV- I - l))x + (1 - pf = 0. (46) 

The physically reasonable solution of Eq. (146]) is x < 1, that corresponds to decreasing 
probability of finding the particle with the distance from the last burnt bridge. Then 
we obtain 

x = 1 +p(N - 1) + ^p\l -l)(N-l-l)- -VS (47) 
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with 



S 



2 + 2p(N- 1) + p 2 {l- l)(N-l- 1) -4(1 -p) 2 . 



(48) 



According to Eq. ( T4T1) . x lies within the unit interval [0, 1] with x = for p = 1 and 
x = 1 for p = 0. We observe that Eq. (146]) and its solution Eq. (1471) are invariant under 
the change I —> N — I, which implies that Eq. (T4"TI) works not only for functions -Rj [Eqs. 
(1391) . (HOI) ], but also for functions Uj [Eqs. (H5j)]. It should also be pointed out 

that Eq. (T46l) cannot be obtained from the recurrent formula (134"]) as it leads to a trivial 
identity. 

In order to find the function Rq from the renormalization condition (132]) . we need 
to eliminate the function Uq. This is done with the help of Eq. ( J33l) [or alternatively 
using Eq. (1361) ]. reflecting the fact that recurrent relations for Rj and those for [/j are 
coupled through Si and S 2 . Comparing Eq. (133]) with (139]) for fc = 0, it follows that 



B{0)=pE 1 



P 



RkN-i + X! UkN-i-i 

k=l k=l 



(49) 



While (1491) can be used directly to express [/q in terms of Rq, it is convenient to combine 
it with another equation, 



5(0) =p 



RkN-i + RkN+i-i 

.fe=l fc=0 



(50) 



to obtain a simpler expression connecting ?7 an d -Ro- To prove the equation (|50|) . one 
has to express RkN-i and RkN+i-i with the help of Eqs. (1401 and (|39l) . correspondingly, 
and with the function B(k) given by ( 1411) . This leads to proving the following expression, 

pR 



RnC 



1 — X 

where we define a new parameter C 
B(k) = Roe ak C, 

which can be written using Eq. (f4"Tl) . 

1 - x + p(N - I) 



{l-lC+(N-l-l)\p + C(-l+p-pl)] + l-(l-l)C},(51) 



C 



(52) 



(53) 



i + (n -i)[i-p+ P iy 

It can be shown that Eq. (15T!) is equivalent to equation (|46|) for x, and therefore (151!) is 
correct, and so is (1501) . Then comparing Eqs. ( 1491) and ( 1501) it follows that 

oo oo 

^ UkN-l-l = RkN+l-1- (54) 
fe=l fc=0 

Expressing U^n-i-i and RkN+i-i according to Eqs. (T42l) and (1391) . we obtain 

U kN -i-! = Uox^l -(N-l- 1)C], (55) 

where we defined another parameter C, 

B(k) = U e ak C, (56) 
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which is equal to [see Eq. 



C = l-x+pl 

N-l + l[l-p + p(N-l)]' { ' 

Combining these results with 

RkN+i-i = RoAl ~ (l ~ 1)C] (58) 
and summing over the index k in (154)) yields 

Uo = i-Q-W „ (59) 

l-(N-l-l)C V ' 

Eq. ( 1591) is the relation between Uq and Ro that is needed in order to determine the 
function Rq. It should be noted that the same result can be obtained by using Eq. (|36|) . 
i.e. B(0) =pL 2 . 

Now we can find Rq from the normalization condition (132j) . which can be rewritten 

as 

ri-i jv-«-i jv-j-i j-i 

1. (60) 



E 

fc=0 



E RkN+i + ^ RkN+l+i + UkN+i + ^ UkN+N-l+i 

.i=0 i=0 i=0 i=0 

We also rewrite functions RkN+i, RkN+i+i and UkN+i, UkN+N-i+i according to Eqs. f|39|) . 
(HOj) . (142]) and (1151) . with functions and JB(A;) given by expressions flU) and flSJ), 

i2 fcJV+i =i?ox fc [l-iC]; (61) 

RkN+i+i =R x k {l-lC + i{p + C(-l+p-pl)}}; (62) 

CW = [W[i -*<?]; (63) 

UkN+N-w = U x k [l — (N — l)C + i{p + <7(-l + p - p[N - /])}], (64) 

where C and C are provided by and (I5"7j) . We note that Eqs. (1631) and (1541) 
correspond to Eqs. (ISTl) and (162]) with I —> N — l change (since the same correspondence 
exists between C and C). The normalization (IBT)]) implies that 

//n f /--/(/-l)C + (iV-/)(l-ZC) 



1 - x I 2 
+ 1 (iV - l)( N - I - l)[p + C(-l + p - pi)]} 

+-^_ - / - I(iV - /)(jv - I - \)C + /(l — [N — l)C) 
1 — x I 2 

+±J(J - l)[p + (7(-l + p - p[JV - /])] } = 1 (65) 
and, introducing parameter B such that Uq = BRq, we obtain from Eq. (159|) 

i» = '-('-Dg., (66) 

l-(N-l-l)C 
Then from (I65p the expression for Rq can be derived, 

#0 = 2(1 - x)/ {21 - 1(1 - 1)C + 2(N - 0(1 - IC) 

+(N-l)(N-l-l)\p + C(-l + p - pi)] 

+5 {2(iV — I) — (N — l)(N — I — 1)C + 2/(1 — (N — l)C) 

+1(1 - l)[p + +p - p(iV - /))]}} . (67) 
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To find the mean velocity of the dimer molecular motor we utilize the formula 
discussed in [9], 

oo oo 

V = J2( U 3 - W j) R j + ^2(U-m- W m )U m . (68) 

j=0 m=0 

The physical meaning of this expression is the fact that the velocity can be viewed as a 
sum of the terms for each possible state of the system as shown in Fig. 8. Due to the 
cancellation of forward and backward rates on all sites except the origin, we obtain the 
expression, 

V = R + U = (1 + B)Ro, (69) 

with functions R and B given by (JBTjl and fl66l) . Eq. fl69l provides the exact formula 
for the dimer's velocity V(l, N, p) (or V(l,c, p)) in the continuous-time forward BBM. 

It can be shown that the derived expression for the velocity agrees with previously 
obtained exact formula in the special cases. First, for p = 1 the parameters needed to 
calculate V are the following, 

. = 0, C=i, e=j±- v B = !Lzl, (70) 
and Eq. (169]) yields 

v ^ = 1) = ,(, + i ) + { N N i)(N-, + iy < 71) 

which reproduces the result (JTj) obtained for p = 1 by using Derrida's method. Another 
special case is when the bridges are shifted by I = N/2 lattice spacings with < p < 1 
(discussed in the preceding section). We checked numerically for complete ranges of c 
and p (with various fixed p and c values correspondingly) and in this way verified that 
Eq. ( !69l produces the correct results. 

The formula (I69p is derived for non-zero shifts, I = 1,2, ■ ■ ■ , N — 1, as discussed 
above. We now consider the case of I — 0, when a dimer crosses pairs of bridges on the 
upper and lower tracks at once. Instead of trying to apply the kinetic scheme in Fig. 8 to 
this case, we observe that a pair of bridges on parallel lattices crossed simultaneously by 
a dimer is equivalent to a single "effective" bridge on the single track but with a different 
burning probability p e ff- Thus the motion of the dimer along two parallel tracks can 
be viewed as hopping of a new effective particle in single-lattice BBM with effective 
burning probability. The velocity for a single random walker in the continuous-time 
forward BBM was found in [9], 

V(N,p cS )= 2Vc l (72) 

^4(iV-l) 2 + 4p eff iV 

The effective burning probability can be found from the following arguments. Both 
bridges can be burnt simultaneously with the probability of p 2 , or only one bridge can 
be destroyed at each lattice, producing 

p eS = 2p{\ - p) + p 2 = 2p-p\ (73) 
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Figure 9. (a) Ratio of dimer's and monomer's velocities as a function of the 
concentration of bridges c for the fixed value of the shift I for p = 0.4. (b) Ratio of 
dimer's and monomer's velocities as a function of the shift I for the fixed concentration 
of bridges c for p = 0.5. 



and Eq. (1721) leads to the following expression, 

V(l = 0,N,p) = , ^P-P 2 ) (74) 

'{2p - p 2 ) 2 {N - l) 2 + A(2p - p 2 )N 



We note that, as expected, p e s > p for < p < 1, and p e ff = P if P = or 1. Thus 
for < p < 1, the dimer's velocity V(l = 0,N,p) always exceeds that of the monomer 
at the same conditions. It is also interesting to note that Eq. (I74"|) can be obtained by 
direct substituting I = into expression ( 1691) that was obtained for / > 1 cases. 

The result (1691) allows us to compare velocities of the dimer and the monomer 
molecular motors in BBM for the entire range of burning probabilities. The ratio of 
velocities at different conditions is plotted in Figs 9-11. The dimer's velocity V2(7,c, p) 
is calculated from Eq. ( l69l) . while the velocity of the monomer V\(c, p) is given by Eq. 
( TT7j) . As shown in Fig. 9a, the ratio V2/V1 shows a non-monotonic behavior as a function 
of the bridge concentration with a maximum that approaches 2 (from below in the limit 
of very large shift distances as long as the burning probability p is not too small. Also 
the ratio of velocities as a function of the shift distance I (for fixed c) has a maximum 
at Ic = 1/2: see Fig. 9b. This behavior is similar to the case of p = 1: compare Fig. 
3 and Fig. 9. The effect of the burning probability on velocities is illustrated in Fig. 
10. Increasing the burning probability leads to the increase in the ratio of the dimer's 
and the monomer's velocities for the intermediate values of c and I, while for very large 
and very small concentrations, as well as very small and very large shifts, the trend is 
reversed. A special case of / = is illustrated in Fig. 11. In this case the maximum 
value of V2/V1 is reached at c = 1 (when every site is a potential bridge), and for all 
concentrations V%jV\ increases as p decreases. The ratio V2/V1 varies from 1 for p = 1 
to \pl in the p — > limit (provided that p<c). 
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Figure 10. (a) Ratio of dimer's and monomer's velocities as a function of the 
concentration of bridges c for the fixed value of the shift I = 5 at different burning 
probabilities, (b) Ratio of dimer's and monomer's velocities as a function of the shift 
I for the fixed concentration of bridges c = 0.1 at different burning probabilities. 
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Figure 11. Ratio of dimer's and monomer's velocities as a function of the 
concentration of bridges c for the fixed value of the shift I = at different burning 
probabilities. 



4. Discussions 

Our theoretical calculations show that the rigid dimer molecular motor always moves 
faster than the monomer in continuous-time BBM. However, the ratio of velocities never 
exceeds 2 (see Figs. 3, 6a, 9, 10 and 11). The behavior of the diffusion constant is more 
complex, as indicated in Figs. 4 and 6b, although the ratio D2/D1 is always between 0.9 
and 1.31. It is interesting to note that for some sets of parameters fluctuations of the 
dimer are smaller than fluctuations of the monomer particle at the same conditions. This 
observation can be understood using the following arguments. The fluctuations of the 
particle in BBM are reduced near the already burned bridges because the particle cannot 
cross them. In the case of the dimer the effective concentration of weak links is larger 
than for the monomer, and the dimer molecular motor more frequently encounters hard 
wall from the left, leading to decrease in the dispersion for some range of parameters. 

More interesting is the behavior of effective forces exerted by the dimer particles 
in BBM presented in Figs. 5 and 7. The dimer always generates larger force than 
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the monomer particle, but in most cases the ratio F 2 /Fi is less than 2. The situation 
is different for the special case of I = 1 and c = 1/2 with the burning probability p 
close or equal to 1. Here the effective force exerted by the dimer becomes very large. 
This behavior can be easily understood for p = 1 burnt-bridge models. In this case the 
motion of the dimer on two parallel lattices can be viewed as the motion of the effective 
particle on a single lattice with larger concentration of bridges - as shown in Fig. 2. 
Because N = 1/c = 2 and I = 1 every link in the effective single lattice is a potential 
bridge. Then the random walker always move forward, and the backward transition 
rate is zero, leading to infinite effective exerted forces for the dimer molecular motor 
[see Eq.flEl])]. 

Computer simulations of the dynamics of clusters of collagenase proteins [6J predict 
that a rigid dimer can generate large forces. In collagen fibers the cleavage sites are 
separated by the distance A = 300 nm, and fibers are shifted by approximately 60 nm 
from each other [HIS]. Although the step-size for the collagen motor protein is unknown, 
one can reasonably take it as ~ 10 nm, the value typical for most known motor proteins 
and consistent with the sizes of these enzymes [U 0]. Then it would correspond to 
iV ~ 30 (or c ~ 1/30), and I ~ 6. Experimental and theoretical estimates of the burning 
probability are between 10 and 30 % [SI [9] • However, even if one assumes that p = 1 
at these conditions our analytical calculations [from Eq. (fl5l) ] produce F 2 — 0.045 pN, 
which is still significantly smaller than predicted from Monte Carlo computer simulations 
for dimers [6]. Thus our analysis suggests that the dimer collagenases, in contrast with 
earlier claims [6], probably cannot produce stall forces comparable to ATP-driven motor 
proteins. 

5. Summary and conclusions 

The dynamics of rigid dimer particles moving along parallel molecular tracks in BBM 
is investigated theoretically using several analytical methods. When the probability of 
burning is equal to 1, we utilized Derrida's method [15] for calculating explicitly dynamic 
properties of the system, and it allowed us to compute the effective forces exerted by the 
dimer molecular motor. For p < 1 the situation when the bridges on the parallel tracks 
are shifted by half of the period, i.e., I = AT/2, is also investigated in detail by mapping 
the system with two tracks into the motion along the single lattice with the distance 
between bridges being N/2. More generally, for p < 1 and arbitrary shifts between the 
tracks, our analysis produced exact analytical expressions for the velocity of the dimer. 

The dynamic properties and the force generated by the dimer are compared with 
corresponding properties for the monomer molecular motor. In all cases it was found 
that the ratio of the dimer's and the monomer's velocities can never exceed 2 (which is 
the limit in case of small concentration of bridges c or large values of the shift /). We 
found that the force generated by the dimer is substantially larger than that generated 
by the monomer only for a special case of c = 1/2 and 1 = 1. For other sets of parameters 
the force exerted by the dimer is at most only slightly more than twice as large as that 
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exerted by the monomer. This occurs for both p = 1 and p < 1 (with I = N/2), 
in the latter case the significant difference between the forces generated by the dimer 
and the monomer at c = 1/2 is present only if the burning probability p is close to 
1. This result is in contradiction with earlier theoretical claims based on Monte Carlo 
computer simulations [6] . These theoretical predictions are explained by the increase in 
the effective concentration of bridges. 

Although our theoretical approach provides exact and explicit expressions for 
dynamic properties of dimer molecular motors in BBM, there are several features of 
the system that should be included in order to obtain more realistic description of 
dynamics of motor proteins that interact with their tracks. It was assumed that two 
motor proteins are tightly bound in the dimer, while in the cell the enzymes must have 
some flexibility. In addition, the interaction between motor proteins in the cluster was 
neglected, although it might lead to the increased dynamic efficiency, as was discussed 
earlier for different motor protein systems [13]. Also, it is not clear what effect on the 
dynamics of molecular motors will have random or non-uniform distributions of bridges. 
It is important to investigate further these features theoretically and experimentally to 
uncover the fundamental mechanisms of motor protein transport. 
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